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1 Introduction 


Regularization by dimensional reduction, a modification of the conventional 
dimensional regularization, was proposed [lj] as a candidate for an invariant 
regularization in supersymmetric theories. The idea is very simple. Super- 
symmetry is only possible in certain integer dimensions (four, to be specific). 
Therefore, the vector and spinor algebra in the numerator of Feynman di¬ 
agrams should be retained four-dimensional. On the other hand, regular¬ 
ization of momentum integrals is achieved in TV = 4 — 2e dimensions. The 
possibility of canceling the squared momenta in the numerator and the de¬ 
nominator, crucial for maintaining the gauge invariance, requires then that 
the momenta form an TV-dimensional subspace in four dimensions. 

Such a dimensional reduction from a finite-dimensional space into the 
formal space of the regularization proves to be mathematically inconsistent 
[§]. The regularized space should be recognized as infinite-dimensional be¬ 
cause antisymmetrization over an arbitrary number of indices does not give 
the identical zero. A consistent modification || disables any strictly four¬ 
dimensional objects, like the totally antisymmetric e tensor, thus distorting 
the algebra above a certain order. On the other hand, the original version, 
in spite of being globally inconsistent, can provide unambiguous results until 
antisymmetrization over five indices actually comes into play (in particular, 
the evaluation of the quark propagator in the quantum chromodynamics may 
became ambiguous not earlier than at the four-loop level). Besides, it is very 
convenient technically, to have the usual four-dimensional algebra of spinor 
matrices (including 75 ) and even their Fierz rearrangement j3J. This is why 
the regularization by dimensional reduction was very soon applied to non- 
supersymmetric calculations || in the Standard Model. It was also used, to 
verify the consistency of the 75 prescription in a three-loop QCD-electroweak 
calculation [|BJ. 

In the present paper we demonstrate the two-loop calculation of the pole 
masses of quarks in quantum chromodynamics in the framework of the reg¬ 
ularization by dimensional reduction. The corresponding calculation in the 
conventional dimensional regularization was done in refs. 0, | 8 j . The leading 
approximation of the present calculation in dimensional reduction (that ig¬ 
nored contributions of heavier quarks and masses of lighter quarks) was used 
in ref. 0 , to verify the 75 prescription. 

Here we evaluate the contributions of the heavier and lighter quarks by 
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means of the asymptotic expansion, which is of a special methodical interest. 
While the structure of the large-mass expansion caused no particular doubts, 
being defined by the universal Euclidean rules 0, the expansion in a small 
mass is more complicated. Its definition needs an on-shell infrared extension 
which we construct explicitly for the present case. 


2 The choice of the renormalization scheme 


Let us discuss the dimensional reduction in some more detail. As the mo¬ 
menta become N -dimensional, 4-vectors naturally split into true iV-vectors 
and so-called e scalars which fall into a complementary subspace of dimen¬ 
sion 4 — N — 2e, orthogonal to the momenta. The £ scalars are nothing 
else but matter fields. Their presence is the only difference from the con¬ 
ventional dimensional regularization. Order by order in perturbation theory, 
contributions of £ scalars are equivalent to finite counterterms, that is, to a 
change of the renormalization scheme. When performing renormalizations, 
one should remember an important fact. Formal 4-covariance (valid at the 
stage of generating the ‘bare’ diagrams) may be broken by counterterms. 
Renormalizations of £ scalars and of their interactions (including the pole 
contributions in £ ) are not identical to those of the vectors. 

Moreover, quantum corrections may also generate a mass for £ scalars. 
There is an arbitrariness in choosing the renormalization scheme for this mass 
0 - A consistent way is to choose the finite £-scalar mass counterterm so that 
the pole (and renormalized) mass of the £ scalars be zero. The £-scalar held 
renormalization is left minimal. Then no additional dimensional parameter 
is ever involved in renormalizations while the polynomial renormalization 
of other masses, independence of renormalizations of dimensionless coupling 
constants on masses, and gauge independence of all /3 functions are retained. 
We follow that scheme in the present calculation. However, this may not 
be the only way. For example, in a softly broken supersymmetric Yang- 
Mills theory it would be natural trying to restore the maximum symmetry 
by relating the mass of £ scalars to the gluino mass as a special solution to 
the renormalization-group equations in the minimal subtraction scheme. 

One should also bear in mind technical complications that are inevitable 
when performing a massive calculation up to a finite part with an additional 
mass on previously massless lines. 


2 




Other renormalizations are done minimally, that is by subtracting only 
poles in e. We use a modified MS definition II dividing each loop by 
(47r) e r(l + e) rather than multiplying by T(1 — e). This definition is more 
convenient in calculations with masses because a simple massive loop is then 
a pure rational function of N without any T functions. It would be proper 
to call the scheme MMS (massive minimal subtractions). 

3 The pole mass 

The singularity mass of a particle is a physically meaningful quantity 0 . We 
restrict ourselves to perturbation theory only and do not analyze the exact 
nature of the two-point function singularity which may involve a branching 
point. In the present paper we accept the conventional term ’pole mass’ 
which refers to the point where the real part of the inverse propagator turns 
into zero. 

Being a physical quantity, the pole mass is renormalization- and gauge- 
invariant. We verify the gauge invariance by performing our calculations in 
the arbitrary covariant gauge. 

The pole mass of a quark mp is defined as a formal solution for p (in the 
Minkowski metric) at which the reciprocal of the connected full propagator 
equals zero: 


p — m — £(p, m) = 0, (1) 

where E(p, m ) = p App 2 , m) +m B(p 2 , m) is the one-particle-irreducible two- 
point function (including the i factor for one of its legs); m may stand for the 
bare or renormalized mass, mp or mp, depending on the prescription used in 
evaluating £. The solution to eq. ([!]) is sought order by order in perturbation 
theory. To two loops 


mp — m + T>i(m, m ) + S 2 (m, m ) + £i(m, m ) E'^m, m ) + 0(S 3 ), (2) 

where is the L-loop contribution to £, and the prime denotes the deriva¬ 
tive with respect to the first argument. 

According to eq. (|2|) , technically, we need to evaluate propagator-type di¬ 
agrams up to two loop on shell. For the one-loop diagrams, also the derivative 
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with respect to the momentum is necessary. It can be expressed through the 
derivative with respect to the mass because the total dimension of the dia¬ 
grams is known. Most of these calculations can be performed with the aid of 
the SHELL2 package |L2| intended to analytically compute propagator-type 
Feynman integrals that involve a continuous line of one mass and the external 
momentum on the mass shell of that mass up to two loops. We use a modified 
implementation of the package, advanced in the scope and efficiency. The 
original algorithm [12|] for reducing scalar numerators is only applied to the 
first powers of the numerators. In other cases we use the strategy described 
in ref. [O] of employing a recurrence relation for a neighbor triangle. 



Figure 2: Further two-loop diagrams, contributing to the pole mass of a 
quark, which involve no other nonzero masses. Solid lines correspond to 
quarks, wavy lines to gluons, and dashed lines to the Faddeev-Popov ghosts. 


At the two-loop level there are 6 diagrams contributing to the pole mass 
of a quark. They are presented in fig. |l](a) and fig. §. 

In individual diagrams there may be on-shell infrared divergences. How¬ 
ever, the complete pole mass is infrared-finite [j7|. It only involves ultraviolet 
poles in e, which are removed after introducing ultraviolet counterterms. 
The counterterm can be evaluated in some other infrared-secure off-shell 
configuration of momenta and masses. The renormalization of the Yukawa 
charge for e scalars and their (nonminimal) mass counterterm are calculated 
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separately. The point is that in a nonsnpersymmetric theory, like QCD, 
the (ultraviolet divergent) renormalization constants for the Green functions 
of e scalars are quite different from those of the vectors (the latter being 
restricted by gauge invariance). Renormalizing a physical quantity, as the 
pole mass, we can accumulate the difference in the bare charge counterterms 
alone, ignoring the field renormalizations, if we evaluate the contributions of 
counterterms for the whole sum of the diagrams. 

4 The rules of asymptotic expansions 

Let us concentrate on the diagram with a quark loop [fig. [T](a)] . When the 
mass of the quark in the loop mi is different from the mass of the external 
quark m, this diagram cannot be calculated by the SHELL2 algorithms. In 
the conventional dimensional regularization, the result involves complicated 
transcendental functions We are going to demonstrate how to evaluate 
the diagram by means of expanding it in the ratio of the masses. If the 
mass in the loop is large as compared to the mass of the propagating quark, 
mi m, then the structure of the asymptotic expansion |J established in 
case of purely Euclidean expansions does not require any modifications. The 
rules are as follows. The expansion is a sum over ’ultraviolet’ subgraphs 
of the diagram. An ultraviolet subgraph must contain all lines with large 
masses, the points were large external momenta (if any) flow in/or out. The 
large momenta ought to go only through the ultraviolet subgraph and obey 
the momentum conservation law. And last, an ultraviolet subgraph should 
be one-particle irreducible with respect to lines with small and zero masses 
although may consist of several disconnected parts. An ultraviolet subgraph 
is Taylor-expanded in its small parameters (external momenta and internal 
masses) and then shrunk to a point and inserted in the numerator of the 
remaining Feynman integral. 

The large-mass expansion of the diagram, fig. [|(a), involves only two ul¬ 
traviolet subgraphs: the whole graph and the loop of two lines with the heavy 
mass. As result of the expansion two-loop bubble integrals and propagator- 
type integrals on shell arise. 

The case mi C m is more difficult. Naively applying the standard ex¬ 
pansion, we would encounter actual infinity for the subgraph that includes 
only the internal ‘on-shell’ line with the mass m. Thus, the expansion of this 
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ultraviolet subgraph needs some extended definition. 


5 The scalar prototype 

It is convenient to perform further considerations in terms of the scalar pro¬ 
totype of the diagram [fig. |l](b)] in Euclidean space (after the Wick rotation). 
The prototype involves arbitrary integer powers of the scalar denominators 
Cl = k\ + m\ on the lines. Their powers j l are called indices of the lines. 
The mass-shell condition for the external momentum now is p 2 = — m 2 . Any 
scalar products of the momenta in the numerator are reduced to changed 
powers of the scalar invariants in the denominator. Thus, the indices may 
sometimes become negative. In particular, an auxiliary line, of the prototype, 
line # 3, is generated with an always non-positive index j 3 < 0. From the 
viewpoint of the R operation or asymptotic expansion, line 3 represents just a 
vertex (local in co-ordinate space). It is convenient to choose = m 2 + mf, 
so that pk\ = |(ci — C3). Other masses are m 2 = 0, m 4 = m, = m±. 

The index of the auxiliary line j 3 can always be reduced to zero by 
means of recurrence relations. The recurrence relations are derived by the 
Chetyrkin-Tkachov method of integration by parts [[b|] . We use a shorthand 
notation { XYZ } of ref. [[13] to denote the relation for the triangle formed of 
lines # A", Y, and Z (the latter is optional, a degenerate case being allowed): 


r d N k x 

/ JxJyJz 
C X C Y C Z 


(n -2jx-j Y -jz+jx 


2m\ 

cx 


+ jv 


rn 2 x + my — 'm X Y 


cy 


+ C X Y - C X 


+jz 


m x + m 2 z - m 2 xz + Cxz ~ c x 
c z 


= 0, 


{XYZ} 


where a double index XY refers to the line that starts at the point where 
lines A" and Y meet. For an external line on the mass shell, the value of Cl 
is equal to zero. Expressing some term of eq. {XYZ} through other terms, 
we obtain a relation between Feynman integrals of the same prototype but 
with varied indices of the lines. 

In order to eliminate a numerator j 3 <0, we apply the following relations. 
If J 5 7 ^ 1, we solve equation eq. {425} with respect to the c 3 /c 3 term. Other 
wise, if J 4 7 ^ 1, the C 3 /C 4 term of eq. {524} is used. If j\ 7 ^ 1, the C 3 
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numerator is reduced using eq. {315} + {425}. In case j\ = j 4 = js = 1 we 
apply eq. {135} + {245}, solved with respect to the free term, to create a 
denominator for which the above relations are applicable. 

At j ',3 = — 1 , one can derive a simpler formula which is more efficient than 
the above algorithm. First of all, C 3 = c 4 — 2 pk\. On integration over k\ 
with just lines 1 and 5 involved, by Lorentz covariance, a single component 
of ki^ turns into k 2fl with the coefficient kih 2 /k 2 . Hence, pki is equivalent 
to {pk 2 ){kik 2 ) / k 2 = j(c 2 — c 4 )(ci + c 2 — c 5 )/c 2 . Thus, the first power of the 
numerator on line 3 can be replace as c 3 = {[c 4 + c 4 + c 5 + (ci — c 5 )c 4 /c 2 — c 2 ]. 

6 The small mass expansion 

We assume that the general i?-operation-like structure of asymptotic expan¬ 
sion is retained in the case of the on-shell expansion in a small mass. The 
set of ultraviolet subgraphs of our diagram (after eliminating line 3) is: {1, 
2, 4, 5}, {1, 2, 4}, {2, 4, 5}, {4}. The first three subgraphs do not cause 
any technical complications in their evaluation. The whole graph becomes 
a two-loop on-shell propagator-type integral which is a particular case of a 
SHELL2 prototype with two massive lines, one of the latter being contracted. 
The two succeeding subgraphs both yield a product of a one-loop bubble in¬ 
tegral by an on-shell propagator type integral. The last subgraph {4} is 
infrared-singular, its formal Taylor-expansion being infinite on mass shell. 
Let us try to derive an extended definition for the asymptotic expansion of 
this subgraph. We need to expand the propagator l/c 4 = 1 /(—2 pk 2 + k 2 ) in 
the small momentum k 2 . It is reasonable to make use of the expansion 


l/(—2 pk 2 + kl) = Y.(kl) n K-ipk 2 + 0)" +1 , 


( 3 ) 


n=0 


where +0 keeps the correct causal (Euclidean) rule of going around the sin¬ 
gularity at zero which appears as we perform the on-shell expansion. 

The even powers of the scalar product in eq. (|3|), for which +0 is inessential 
in the context of dimensional regularization, can be worked up by the formula 


/ 


d N h (pk 2 ) 2n f(k 2 2 ) 


r(7V/2)r(l/2 + n) 
T(7V/2 + n)T(l/2) 


/ 


d N k 2 e v 2 ki) n f{ki ), (4) 
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where / is an arbitrary function which depends on k 2 through k 2 only. The 
formula can be proved by induction based on the Lorenz structure of the 
result, for positive integer values of n and then analytically continued to 
arbitrary complex values. It is regular for negative integer n. We assume 
that j:i has been brought to zero by recurrence relations before the expansion. 
Then, subgraph {1, 5} after integration over k\ depends only on k? 2 . This fact 
permits us to apply eq. (f|) to even negative powers of pk 2 in eq. (Q). The 
integral can then be explicitly evaluated and transformed to the form 

r d N k\ d N k- 2 

J % N T 2 (3 - iV/ 2 )ciV 2 V 5 5 (pA; 2 )^ “ 

(mj) N -h-j2-h-j T ( 1 / 2 )T(N/2 - j 2 - j) 

(-p 2 ) i r(j 1 )r(j 5 )r(i /2 + j)r(N/2 - j) 

T(ii + h + j - N/2)T(j 5 + j 2 +J- N/2)T(j 1 + j 2 + j 5 + j - N) 
T(ji + is + 2 j 2 + 2 j — N ) 

Odd negative powers reveal on-shell infrared singularities of the propaga¬ 
tor, which ought to generate some counterterms as a part of the asymptotic 
expansion. They become manifest as we remember that 

1 / 0^2 + 0) 2n+1 = 1 /(pk 2 ) 2n+1 + 5 ( 2 n) (pk 2 ) ( 6 ) 

( 2 n)! 

for purely imaginary p as we have it. The first term on the right-hand side 
of eq. ([]), which generates the integral in the sense of the principal value, 
yields just zero with a function of k 2 . The second term picks out exactly the 
2n’th coefficient in the Taylor expansion of the rest of the diagram in pk 2l 
that is the component of k 2 parallel to p. This is the kind of infrared coun¬ 
terterms, only partially local in momentum space, which are characteristic of 
the non-Euclidean asymptotic expansion. The delta function itself takes off 
one integration over the corresponding component of k 2 , leaving the integral 
in dimension N — 1. 

The seemingly independent treatment of the odd powers of the scalar 
product in the denominator by eq. @ is in fact firmly bound to formula (Q) 
for the even powers in the general framework of dimensional regularization. 
The point is that the latter always implies an auxiliary intermediate analytic 






regularization in the line indices []. One seeks a correlated region of all the 
parameters, the powers and the space-time dimension, in which the integral 
under consideration is well defined and analytic in the parameters. Then the 
resulting function is analytically continued to the whole complex plane, and 
eventually, the limit is taken in the indices to the values that we started with. 
As a rule, the regular limit does exist while the space-time dimension is kept 
noninteger. According to this schedule, we have evaluated the integral on 
the right-hand side of eq. (££) for the relevant diagram, using the power of the 
scalar product as an additional regulator rather than a natural number and 
taking the limit n —> — j, eq. Having identically transformed the result, 
we can try to take the limit for half-integer values of j as well. The limit does 
always exist and exactly coincides with the value that we get from / d N ~ 1 k 2 
according to the eq. (|). It is interesting to note that the odd coefficients of 
the expansion turn out to be ultraviolet-finite, that is allow a smooth N —> 4 
limit. Moreover, the limiting odd series happens to terminate because of a T 
function in the denominator. 
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Figure 3: The one-loop diagram where the asymptotic on-shell expansion 
in m 2 -C 777.4, V 2 = ~ m \i is performed. 

The use of eq. ([6|) gives us an alternative way of achieving the same result, 
without specifying the weight function and avoiding the necessity to evaluate 
the indeterminacy that arises at half-integer negative n in eq.([|) because of 
a singular coefficient. We have verified both procedures in yet another case 
when the integration in eq. (jj) can be performed analytically in terms of T 
functions. This is the scalar one-loop diagram with one light mass and a line 
on mass shell (fig. ]3|) Q. For odd expansion coefficients we obtain the same 

3 L.V. Avdeev is grateful to Prof. K.G. Chetyrkin for a clarifying discussion on the 
subject in 1983. 

4 On preparing the present paper, we became aware that a general method of con¬ 
structing typically Minkowskian and on-shell asymptotic expansions was proposed by V. 
A. Smirnov |l5| and demonstrated by one-loop examples (including (fig. )) and by a 
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results either by analytically continuing eq.(f|) or by reducing the dimension 
of integration to N — 1 via eq. (|). The complete expansion, which includes 
the contributions of subgraphs {4} and {2, 4}, perfectly reproduces the exact 
integral in the Feynman parameter 


(/I 


2\e 


d N k 


7 r v/ 2 r (3 _ N/2)(k 2 + m 2 2 ){-2pk + k 2 ) 


2 - In 


= - + 2 - In 

£ 


~P 

, P? 

'_p2 

, k 2 


— r 2 lnr — rv/4/r 2 — 1 arctan J A/r 2 — 11 + 


/ r 
— 7r r + — + 


128 


— r 2 (ln r — 1 )-h 

V ; 12 


0{e) 

(7) 


where r 2 = m 2 /(—p 2 ). It is worth pointing out the fact that analytic contin¬ 
uation recovers just the causal rule of getting around the singularity, as in 
eq. ( 0 ), which is automatically implied in the Euclidean space where k 2 > 0 
after the Wick rotation. 

Relations (|3j)-(b) are sufficient [with a trivial generalization of eq. (|3) by 
binomial coefficients as index j 4 is arbitrary] for constructing the asymptotic 
expansion terms for the ultraviolet subgraph {4} of the prototype, fig. |l](b). 
The procedure does not require any modifications for a higher-loop analog 
of fig. |(a) with an arbitrary insertion in line 2 . 

The expansion of ultraviolet subgraph {4} suffers from on-shell infrared 
singularities. Therefore, the expansion is not perfectly local in coordinate 
space (or polynomial in momentum-space). Along with regular Taylor-like 
terms, (& 2 ) n , partially local on-shell infrared counterterms S^ 2n \pk 2 ) are pro¬ 
duced in odd coefficients. At j 4 > 1 , also nonlocal negative powers of k 2 
appear in even coefficients. However, the general structure of the expansion 
as a sum over ultraviolet subgraphs remains valid. 

In case of the conventional dimensional regularization the sum over all the 
ultraviolet subgraphs of the diagram (fig. 0) exactly reproduces the expansion 
presented in ref. ||. 


non-trivial two-loop master-integral j[(| . 
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7 Techniques and results of the calculations 

The calculations in dimensional reduction do not generate any new proto¬ 
types of Feynman diagrams. Just the coefficients before the same scalar in¬ 
tegrals are slightly changed for the bare diagrams. If we use iV-dimensional 
indices, the full contribution of the vectors and e scalars to a diagram with¬ 
out external gluon lines is obtained by replacing N in the numerator of a 
scalarized expression with 4. The contribution that involves £ scalars can be 
extracted by replacing N n with 4” — (4 — 2e) n . 

If we need to separate the proper-vector and e-scalar propagators we 
should proceed as follows. The gluon propagator diagram is scalarized by 
contracting it with the metric tensor either in the formal four-dimensional 
space that contains the (4 —2e)-dimensional vector and e-scalar subspaces, or 
just in 4 — 2e dimensions. According to the Ward-Slavnov-Taylor identities, 
the vector part (for the sum of all diagrams in the relevant order of perturba¬ 
tion theory) ought to be transverse with respect to the external momentum, 
while the e-scalar part is always diagonal as there are no momenta in this 
subspace. Again, we perform the algebra using the formal N-dimensional 
indices. As the external indices are yet free, we should set IV = 4, in order 
to sum up the vector and e-scalar contributions in the internal lines of the 
diagram. Then the external indices are contracted, which produces N n with 
n — 0,1. The e-scalar propagator is then obtained by replacing N n with 
[4 n — (4 — 2e) n ]/(2e). The scalar coefficient before the transverse projection 
operator in the vector propagator can be evaluated by replacing N n with 
(4 — 2e) n /(3 — 2e). The momenta are always assumed to belong to 4 — 2e 
dimensions. 

According to eq. (|2|), the pole mass acquires additive contributions from 
the two-loop diagrams of the quark propagator. The evaluation of the e- 
scalar contribution to the one-loop diagram and its derivative on shell is an 
easy task performed analogously to the corresponding calculation in the con¬ 
ventional dimensional regularization [7], Q. The construction of the small 
mass asymptotic expansion for the diagram of fig. |T|(a) was discussed above. 
Other two-loop diagrams shown in fig. |2] are, in principle, computed straight¬ 
forwardly on the same lines as in the conventional dimensional regulariza¬ 
tion |8|]. The unrenormalized results for individual bare two-loop diagrams 
are summarized in the Appendix. For comparison we also present the corre¬ 
sponding contributions in the conventional dimensional regularization. 
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Here are our results in the framework of the dimensional reduction MMS 
scheme with the zero pole mass for e scalars. We have evaluated the one-loop 
vertex and held renormalization constants: 


Zffv 

Z ffs 


h r,£- 1 

1 - - ^ 


£ L 


(~4 -hi )C a +t;Cf + 0(h 2 ), 

1-—-— C a + (£ + 1) Cf + 0(h 2 ), 


1 - 


1 - - 
£ 1 


h r 


(3 -OC a -N f + 0(h ). 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 


Indices /, v, s denote the type of the particles, fermions, vectors, and e 
scalars. The expansion parameter is h = g 2 /(167r 2 ) = a s /(47r); £ is the 
parameter of the covariant gauge, £ = 1 corresponds to the Feynman gauge, 
£ = 0 to the Landau gauge. We assume an arbitrary simple Lie gauge group 
of colors and an arbitrary irreducible representation for the quarks. In case of 
SU(n) and the fundamental representation, the group coefficients would be: 
C a = n, Cf = (n 2 — l)/(2n). The trace normalization was chosen as Tf = 1/2. 
The number of quark flavors is Nf. The held renormalization constants are 
defined for the one-particlc-irreducible two-point diagrams, including one leg. 
They enter the charge renormalization in positive powers: 


^ = Zf fv Z 2 f Z v = l + ^(^Ec a + ^N f )+0(h 2 ), (12) 

Z y = Z 2 fs Z 2 Z s = l + , l(-2C a -2Cf + N f ) +0(h 2 ). (13) 

As expected, the gauge and Yukawa charges are renormalized differently, 
eqs. (|1^) and ([L3]), because there is no supersymmetry to protect the tree- 
level coincidence of the e-scalar and vector coupling constants. We every¬ 
where replace the Yukawa charge with the corresponding function of the 
gauge charge, that is, exploit a special solution of the renormalization-group 
equations. 

The nonminimal mass counterterm for the e scalars 
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+ 0{e) + Oih 2 ) 


(14) 


m 2 f \ 

7?) 

is necessary to insure that their pole mass is zero. 

The fermion held and mass renormalization constants were calculated up 
to two loops: 


A m 2 = 



2 In 


Z f = 1 


h 


Rev) 


Cj} + 0(h 3 ), 


4 e 

i , 7 + ini + «- i) 2 i 




7 — 

m 


L2e 2 4e 
1 -±3C,-#[(-£ + ”.)C.C, 

+ (-i-sK + (?-^) c R +0(fc3) ' 


CaC, 

(15) 


(16) 


In spite of the fact that there are several masses in the theory, there exist a 
common mass renormalization constant for all flavors. Let us remind that in 
the conventional dimensional regularization the mixing of the masses through 
the diagram of fig. |l](a) was prevented by the gauge invariance which forbids 
a mass term in the vector propagator. In case of the £ scalars, there is no 
such ban in general. But subtracting the e-scalar propagator at zero, we, at 
the same time, get rid of the m 2 /e term in the incomplete R operation. This 
is related to the fact that the subtraction at zero insures the absence of an 
infrared singularity in the renormalized subgraph {1,2,5}. 

The final result for the pole mass of a quark with a given flavor / is 
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v m ' f 


= 1 + h 


771 f 

3 In —f 
M 2 


C f 


* a {[ ; 


1093 11, 2 m 2 f 179, ml 

_I_]jp __In _ L 


24 


+ 


/' 


6 


In ■ 




- 6C(3)-8C(2) + 24C(2)ln2 


C a C f 


+ 


59 9 


,2™) 


m, 


-— + - In- —f + - In —£ 


I 1 


p 


12C(3) + 30C(2)-48C(2) In 2 


Cl 


, 2 m / 13 , m / 
g[-“'V + T ln 1 ? 


2 "" f , + + ( 17 ) 


The sum over flavors represents the renormalized contribution of the diagram 
with the quark loop, fig. [I](a). We know this contribution in three kinds of 
expansions. The small-mass expansion, the construction and meaning of 
which was discussed above, gives 


37 

E(r) = -4C(2) + 2 tt 2 r - 12 r 2 + 2 vr 2 r 3 

6 


+ 


151 . 13, 2 ,2 2 ' 

-4£(2) H-ln r 2 — In 2 r 2 

18 3 

r 1389 9 


1 A 

[-76 

8 , 2 l 

r 4 + 


— In r 

-i 

-75 ~ 

15 


— — In r~ 


-9800 70 

+ C>(r 12 lnr 2 ). 


+ 


r 3988 16 


L99225 315 


In r‘ 
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(18) 


Numerically the most essential correction due to the mass of a light quark 
corresponds to the linear term in the mass ratio. 

The large-mass expansion, performed according to the Euclidean rules, 
yields 
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w N 20 13 2 , 2 2 /76 8 2 \ , 

A(r) =- 1 -In r 2 + In' r 2 — (- 1 -In r 2 )/r 

K " 3 Wr IS )' 


'75 15 

/1389 9 2 \ / 4 / 3988 16 2 \ / 

-(- 1 -lnr 2 j/r 4 - f- 1 -Inr 2 h 

Vqsnn 7n A Vqqook sir )' 


1229 


178408 ' 198 
+(9(r -12 lnr 2 ). 


In r 2 J /r 


'99225 315 

/ 184452 


125050025 5005 


72 , 2 \ 

+ In r J /i 


10 


(19) 


In most calculations, people simply ignore the quarks that are heavier than 
the characteristic energy scale of the problem. They use thus an effective low- 
energy theory with the reduced number of particles. The running charges 
and masses of this effective theory are related to those of the full theory by 
the so-called matching conditions that the running parameters are supposed 
to change their running rate at some value of the renormalization scale /i 
of the order of the heavy-particle mass. Of course, the low-energy effective 
theory is no longer valid at that scale. The matching condition just ensures 
the equivalence of the two theories at low energies where the heavy particles 
ought to decouple, up to power corrections, irreproducible in the low-energy 
theory. The coincidence of the pole mass in the two theories provides the 
two-loop matching condition for the running mass. The charge matching is 
done by means of another physical quantity, like the invariant charge. In 
the MMS scheme (as well as MS), on the one-loop level, relevant here, it is 
done at /i = rrij [[H]. One should bear in mind, however, implicit numeric 
ambiguities of any matching condition, similar to the scheme dependence, 
owing to the fact that higher orders of perturbation theory are unavailable. 
The result, thus, depends on the choice of a physical quantity and on the 
detailed prescription of equating it at low energy. When defining the two-loop 
mass matching, the term of order zero in the ratio of the masses is absorbed 
into the matching point. The actual correction starts from the leading power 
m 2 /m 2 . 

The third available expansion is performed in the difference of the squared 
masses. Such an expansion covers the intermediate region between eqs. (|H|) 
and (|1|). In the present case the intermediate expansion is just a regular Tay¬ 
lor expansion involving only the whole diagram as the ultraviolet subgraph. 
A surprising remark is, however, that generating this trivial expansion is the 
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E (r) 



Figure 4: Three expansions for the renormalized diagram of fig. 0(a) in 
the dimensional reduction MMS scheme. The dotted line corresponds to the 
small-mass expansion up to r 6 , the dashed line to the large-mass expansion 
up to 1/r 6 , the solid line to the intermediate expansion up to (r 2 — l) 3 . 


16 






most time-consuming of the three, while the most complicated small-mass 
expansion is the fastest. 
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( 20 ) 


As the masses of different quarks are quite different, higher terms of eq. (] 2 (i|) 
are of little practical use. The value that is important is 77(1) which defines 
the self-contribution of a quark to its pole mass. The intermediate expansion 
just demonstrates the means of recovering the function E{r ) in the whole 
region of the mass ratio. The matching of the three expansions is demon¬ 
strated in fig. |3|. We terminated the expansions at the same relatively low 
order for the sake of clearness. Keeping more terms, we get a wider region 
around r = 1 where the curves practically coincide. 


8 Discussion 

The results of the three expansions flT8|)-(j20j) differ from the corresponding 
expansions in the conventional dimensional regularization MMS scheme only 
by a constant — j in the zeroth order, which confirms that these are just two 
different renormalization schemes. 

All the three expansions reproduce the same dependence of the diagram 
on /i. This confirms the fact that our renormalization scheme, though being 
nonminimal, remains massless. The latter implies that the renormalization- 
group functions do not depend on masses and other dimensional coupling 
constants but polynomially. The nonminimality of our scheme nevertheless 
reveals itself by the fact that the f3 function of the fermion mass cannot be 
extracted from the first-order pole of the renormalization constant (p~ 6 |) . The 
correct way of defining the mass beta function (valid in an arbitrary renor¬ 
malization scheme) is to extract it from the condition of the renormalization 
invariance of a physical quantity, the pole mass, by differentiating eq. ( 0 ) 
with respect to ln/i 2 . The result 
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( 21 ) 


does not coincide with what has been obtained in ref. [ IS 1 in the minimal 


subtraction scheme, although the renormalization constant (p~6|) , the one-loop 
counterterms, (1^)-( |I3|) , and the one-loop contribution to eq. (|T5|) , agree with 
refs. [|1^, [Hj, which confirms that all the calculations are correct. 

The situation is explained quite evidently indeed. In our renormalization 
scheme the bare mass of the fermion ceases to be a renormalization-invariant 
quantity. The point is that the physical mass involves contributions from 
two mutually correlated bare masses - of the quark and of the e scalars. 
Both do depend on //, to ensure that the physical quantity is invariant. The 
true reason for the difference between the f3 function and the 1/e pole in the 
renormalization constant is thus an extra bare parameter as compared to the 
number of the physical renormalized parameters. In fact, the renormalized 
mass of the e scalars becomes one more independent quantity. However, we 
keep it equal to zero, while in the minimal subtraction scheme it would be 
running as well. Total subtraction of loop corrections to the e-scalar mass, as 
we perform here, seems to be practically the only scheme applicable to phys¬ 
ical calculations (including the finite parts) with masses in nonsupersymmet- 
ric theories, if we are going to stick as closely as possible to the original idea 
of dimensional reduction which naturally suggests the zero tree-level mass 
for the e scalars as well as for the gauge vectors. The minimal subtractions, 
being admissible theoretically and convenient for renormalization-group cal¬ 
culations, would involve an extra independent mass parameter, irrelevant to 
actual physics, however, up to a scheme redefinition of other parameters [EH] 


In principle, all renormalization schemes are, of course, equivalent. The 
difference of our (3 function (^TJ) from ref. |18] is just a renormalization- 


scheme difference which can be compensated for by a finite recalculation of 
the parameters. 

Our renormalization scheme can be related by proper recalculation to the 
standard dimensional renormalization (without the e scalars) as well. A re¬ 
calculation can always be done by equating a physical quantity evaluated in 
the two schemes. Our result © allows us to recalculate the renormalized 
mass at the two-loop level. However, we have to take into account the fact 
that the renormalized charges in the two schemes are also different. To re- 
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late them we should equate another renormalization-invariant quantity like 
the one-loop invariant charge. The necessary relation between the renormal¬ 
ized charge in dimensional reduction and in the conventional dimensional 
renormalization can be taken from ref. ||: 


hoRED — {h 1 + ~C a h + 0(h 2 ) | 


DREG 


Then the mass recalculation looks like 


( 22 ) 


?«dred — 1 — hCf + h~Cf(^ 
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12 2 * 4 ■ 


+ 0(h 3 )]} 


DREG 


(23) 


One could also derive the recalculation formulae from comparing the re¬ 
normalization-group functions in the different schemes. However, then one 
would need to calculate one-loop further. For example, the comparison of the 
two-loop (3 function (pl|) with the corresponding result in the conventional 
dimensional renormalization reproduces the one-loop term of eq. (|23|). 

Acknowledgments This paper was supported in part by RFFI grant 
# 96-02-17531, JSPS FSU Project and the Heisenberg-Landau program, and 
on the part of Kalmykov, by Volkswagenstiftung. Kalmykov is grateful to 
the Physics Department of Bielefeld University for warm hospitality. 


APPENDIX A. Contributions of bare two- 

loop diagrams to the pole 
mass of a quark 

We present the unrenormalized contributions, implying a common factor 
h 2 m(/x 2 /m 2 ) 2e , up to discarded 0(e) terms, in the regularization by dimen¬ 
sional reduction and the conventional dimensional regularization. In the 
former case we use the parentheses, while in the latter case the brackets, in 
the notation of the diagrams. 
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fig- 0 a ] = C f[~^ 
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where mi is the mass of the quark in the loop, and E(r ) is defined by expan- 
). Note that the result in dimensional reduction ( A.l ) involves 


sions (|T8|) - 

a mass-mixing term proportional to rn\. This term is, however, canceled lat¬ 
ter by the ^-scalar mass counterterm. Also, the e-scalar held renormalization 
counterterm affects the 1 /e pole of the renormalized diagram, so that in the 
end the difference from the conventional dimensional renormalization is re¬ 
duced to — j. The subsequent diagrams in the two schemes similarly involve 
differences in the lower-order poles and the finite parts: 
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fig. 0 (c) 
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